Nonmonotonic dj.2 



y2 superconducting gap in electron-doped Pro.sgLaCeo.n 
Evidence of coexisting antiferromagnetism and superconductivity? 
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Recent experiments on Pro.89LaCeo.iiCu04 observe an anisotropic spin-correlation gap and a 
nonmonotonic superconducting (SC) gap, which we analyze within the framework oi a, t — t' — t" — 
t'" — f" — U model with a d^2_y2 pairing interaction including a third harmonic contribution. By 
introducing a realistic broadening of the quasiparticle spectrum to reflect small-angle scattering, our 
computations explain the experimental observations, especially the presence of a maximum in the 
leading edge gap in the vicinity of the hot-spots. Our analysis suggests that the material behaves 
like a two-band superconductor with the d-wave third harmonic acting as the interhand pairing gap, 
and that the anti-ferromagnetic (AFM) and SC orders co-exist in a uniform phase. 

PACS numbers: 74.72.Dn 74.20.Rp 74.25.Dw 74.25. Jb 



Recent Raman scattering ll| and ARPES 
experiments y on electron-doped Pro.89LaCeo.iiCu04 
(PLCCO) find a d,j.2_y2 superconducting (DSC) gap, 
which varies non-monotonically along the Fermi surface 
(FS). The gap increases as one moves avi^ay from the 
nodal direction towards the hot-spot region in the 
Brillouin zone (BZ) where it attains a maximum value. 
It then decreases as one approaches the zone boundary 
along the anti-nodal direction. The involvement of the 
hot-spots suggests that the bosonic pairing possesses 
a magnetic origin. [3|, Ij] If we assume that the system 
continues to remain uniform with doping, and recall that 
the residual gap in a half-filled AFM insulator yields 
two FS pockets on electron doping, these observations 
present a striking conundrum: The leading-edge gap is 
the largest in the momentum region of the hot-spots 
where there are no segments of the FS and quasiparticle 
states lie well below the Fermi energy. 

In this article, we discuss a relatively simple route for 
resolving this dilemma and understanding the behavior 
of the superconducting state of PLCCO. In particular, 
we consider a uniformly doped system with co-existing 
AFM and DSC orders, where the quasiparticle spectrum 
is broadened realistically to reflect the effects of small an- 
gle scattering. We thus obtain in the normal state a finite 
spectral weight at the Fermi energy (Ep) throughout the 
BZ, and when a third harmonic term is included in the 
pairing interaction, the computed leading edge gap along 
the FS reproduces the corresponding experimental varia- 
tions remarkably well. The third harmonic contribution 
possesses the proper symmetry to couple the DSC order 
parameter on the two FS pockets to induce a maximum 
in the leading edge gap around the hot-spots even though 
the hot-spots lie outside the momentum region of the two 
FS pockets. Notably, non-monotonic gap variations have 
been obtained in some earlier calculations!^ Q, but we 
are not aware of a previous study involving a realistic 
model of electron doping with two FS pockets. 

Our study gives insight into the issue of well known 
asymmetry between the properties of cuprates with elec- 



tron vs hole doping, which has been a subject of consid- 
erable recent debate (see e.g., Ref. p). There is grow- 
ing evidence that the normal state of the electron doped 
cuprates can be described as an AFM metal up to a quan- 
tum critical point near optimal doping, and that nano- 
scale phase separations or stripe physics so prominent in 
the hole-doped case arc weak or absent in the electron- 
doped systems. 6, 7] Since our analysis based on a single 
phase model in which the AFM and DSC orders co-exist 
is able to reasonably explain non-monotonic gap varia- 
tions, it supports the scenario that the system remains 
uniform with electron doping not only in the normal but 
also the SC state without the intervention of other orders. 
We model PLCCO as a uniformly doped spin density 
wave (SDW) antiferromagnet with d-wave superconduc- 
tivity. At the mean field level, there is long-range AFM 
order, but when fluctuations are included, the Neel tem- 
perature Tn becomes zero (in the absence of interlayer 
coupling), and the mean- field gap turns into a pseudo- 
gap A* with crossover temperature T*. In the SC dop- 
ing range the FS consists of a necklace of two types of 
pockets: an electron-like pocket near the antinodal point 
(tt, 0), and a hole-like pocket near the (7r/2,7r/2) nodal 
point. A pseudogap near the hot-spot separates the two 
pockets. Our one-band Hubbard model Hamiltonian is 
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where cl (cg^) is the electronic creation (destruction) 

operator with momentum k and spin a. The independent 
particle dispersion with respect to the chemical potential 
Ep is given by[i,|g| 

^j: = -2t[c^{a)+Cy{a)]~U'c^{a)cy{a)^2t"[c^i2a) 
+Cy{2a)] - At'"[c^i2a)cy{a) + Cy{2a)c^{a)] 
-'irc,i2a)cyi2a)-EF, (2) 



with Ci{aa) = cos{akia) and a is the lattice constant. 

The effective on-site Hubbard repulsion Uq is taken to 
be doping-dependentjg. The Mott gap UqS arises from 
a finite expectation value of the staggered magnetization 
S at the commensurate ordering wave vector Q = (tt, tt). 
The d-wave superconductivity including first and third 
harmonics is defined byQ, [ifll 

Ar = A,r + A^r = y"v(k,k'){cl c^ - ) 

k Ifc 3fc A^ ^ ' ^\ k',1 -fc',X' 

k' 

= E ^^9^k - E 5.e^^ E 9.k' (^M^-Pa)' (3) 

with i = 1,3 and g^^: = Cx{ia) — Cy{ia). 

Hamiltonian of Eq. 1 is diagonalized straightforwardly 
|lO| . The resulting quasiparticle dispersion consists of 
upper and lower Hubbard bands (UHB and LHB), each 
gapped via the DSC pairing as: 
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FIG. 1: (color online) (a) Computed quasiparticle spectrum 
of PLCCO at 11% doping. Artificially large values of the 
pairing interaction parameters (Ai = 20 meV and A3 = -6 
meV) have been used to highlight the effect of the SC gap, 
which splits the UHB and the LHB into the pair of bands 
E^ (magenta and red) and E^ (cyan and green), respec- 
tively. Spectral weight is represented by the vertical width of 
the bands, (b) FS map in PLCCO at 30 K (normal state) 
obtained by integrating the spectral intensity of Eq. |7|over a 
small energy window of ± 5 meV around Ef- High intensity 
is denoted by red and low by blue. 



where EK ^ (q)^ + [UqSY and ^± = (^g ± f£+Q)/2. 
The SDW magnetization 5* per site is evaluated self- 
consistently at each temperature via the equations. 
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where f{E) — l/(exp {(3E) -I- 1) is the fermi function with 
B — 1/kBT. The expressions for ar, Br, u- and u- are 
same as Eq (9) and Eq (14) respectively of Ref. [lO|. The 
self-consistent SC-gap equations are 
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The one-particle Green function is 
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where we have included broadening due to small-angle 
elastic scattering [llj, [l^ with the associated renor- 

malization factor, Z = 1 + iT o{Lu)sgn{Lu) / y/ lu'^ — A-. 
Tq{uj) is the normal state scattering rate, discussed be- 
low. The corresponding spectral intensity is A(k, lu) = 
-Im[G(k,w)]/7r. 

We fitted the dispersion given by Eq. ^to the ARPES 
resuhs on PLCCO0 at 30 K to obtain the relevant TB 
parameters as follows: t = 0.12 eV, t' = -0.06 eV, t" = 



0.034 eV, t'" = 0.007 eV, & = 0.02 eV, with Ep = -0.082 
eV at 11% doping. [131 The effective Hubbard parameter 
is found to be Uq = 4.15t, which yields a self-consistent 
value of the magnetization S of 0.281 from Eq. [S] 

Fig. 1 clarifies the nature of the quasiparticle spectrum 
and the FS. As Eq. 4 shows, the non- interacting band 
is split via the Hubbard interaction into the upper and 
lower Hubbard bands. As seen from Fig. 1(a), each of 
these bands is further split by the SC interaction to yield 
pairs of bands E^^ (magenta and red) and E^^ (cyan 
and green) which lie symmetrically above and below the 
Ep. [Note that the SC gap has been made artificially 
large in Fig. 1(a) to highlight the effect of the SC split- 
ting.] The spectral weights for various bands (propor- 
tional to the vertical width of shading) are seen to vary 
greatly with k. In effect, the spectral weight of the non- 
interacting band is redistributed between the UHB and 
the LHB through the AFM order, and it is modified fur- 
ther via the SC order at very small energy scales of a few 
meV's. Above the Ep, we see from Fig. 1(a) that most 
of the spectral intensity resides in the UHB branch E^^ 
(magenta) along the (tt, 0) — > (tt, tt) — > i'^/'^i ''^f'^) 
line. At other momenta, the spectral weight generally 
lies below Ep and is carried by the i?^~ band (green). 
Fig. nib) shows a map of the normal state (30K) FS. 
This FS is consistent with the experimental results on 
electron-doped Ndo.87Ceo.i3Cu04 (NCCO) J4, J^j, and 
it can be understood with reference to the band structure 
of Fig. 1(a). The E^~ band (red) gives rise to the (7r,0)- 
centered electron pockets, while the (7r/2, 7r/2)-centered 
hole pocket arises from the E'"^ band (green). 

Fig.|21gives further insight into the nature of the quasi- 
particle spectrum and the associated spectral intensity 
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FIG. 2: (color online) (a)-(c) Quasiparticle spectrum of 
PLCCO with 11% doping in the normal state at 30 K along 
three different lines (see insets) in the (fc^;, fcy)-plane. Spec- 
tral weights are proportional to the widths of lines as in Fig. 
1(a). Blue dots give the corresponding experimental disper- 
sions taken from Ref 2. (d)-(f) Spectral intensity obtained 
from Eq.[7|as a function of binding energy for a series of mo- 
menta where the actual momentum values correspond to the 
experimental points (blue dots) in (a)-(c). 



in PLCCO. The computed energy bands in the normal 
state along the three cuts used in the fitting procedure 
are compared directly with the corresponding experimen- 
tal dispersions in panels (a)-(c). The overall agreement is 
seen to be quite good, although a discrepancy is evident 
in (a) for the anti-nodal cut in that the computed size 
of the electron pocket around (7r,0) is smaller than that 
indicated by the experimental data. [16] Notably, due to 
the presence of the AFM gap, quasiparticle states along 
the hot spot in (b) lie well below the EF^ leading to a 
suppression of the spectral weight up to binding energies 
of about 60 meV. In contrast, along the nodal direction 
in (c), the quasiparticle band crosses Ep^ and as already 
pointed out, gives rise to the (7r/2,7r/2) centered hole 
pockets. 

The spectral intensity computed from the imaginary 
part of the one-particle Green's function of Eq. 7 is shown 
for a series of momenta in the bottom row of Fig. 2. Here 
the cuts in the top and bottom set of panels correspond 
to each other. For example, in (e) the five spectra shown 
refer to the five k-points given by the blue dots in (b) with 
the lowest spectrum in (e) corresponding to the leftmost 
dot in (b). The computations assume a fc-independent 
broadening function of form: ro(a;) — Co[l + {ijj/ujqY], 
which is similar to a form which has been proposed for 
hole-doped cuprates[l2|. Parameter values of Co — 0.1 
eV, ujQ ~ 1.59 eV, with p =3/2 reasonably reproduce 
the experimentally observed broadenings 1^ . We empha- 
size that, for our purposes, the detailed spectral shape 
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FIG. 3: (color online) Spectral intensity computed from Eq. 
|7|in the SC state at 8 K (solid lines) is compared with that 
for the normal state at 30 K (dashed lines) at four representa- 
tive momentum points A-D (indicated by stars in the insets). 
Vertical dotted lines mark the -Ef, which corresponds to the 
energy where the normal state spectral weight falls to half of 
its value at higher energy. The shift (Ashift) of the spectral 
intensity away from the Ef as a leading edge gap opens up 
in the SC state is seen in (a)-(c), while this shift vanishes in 
(d) along the nodal direction due to the d-wave symmetry of 
the gap. 



is not so important. The key is the presence of quasi- 
particle broadening, which allows the development of a 
finite spectral weight at the Ep and the formation of the 
leading edge superconducting gap at all momenta, even 
though the underlying quasiparticle states lie well below 
the Ep at most momenta. 

Fig. 3 compares spectral intensities in the vicinity of 
the Ep in the normal and the SC state at four different 
momenta A-D (see insets). We see that in (a)-(c) the 
midpoint of the leading-edge of the 8 K spectrum (solid 
lines) is shifted by Ashift towards higher binding ener- 
gies with respect to the spectrum at 30 K (dashed lines) , 
while along the nodal direction in (d), this shift vanishes 
due to the node in the d-wave gap. Experiments |2| find 
a leading edge gap Ashift of 2.0 meV along the antin- 
odal direction at the momentum point A (7r,0.8), which 
increases to 2.5 meV at the hotspot B(2.34,1.0); Ashift 
then decreases to 1.6 meV at C(1.85,1.23) and becomes 
zero along the nodal direction at the point D(1.52,1.52). 
The computations of Fig. 3 reproduce this behavior. For 
this purpose, the superconducting parameters are found 
to be: V\ = -76 meV and V3 — 54 meV, leading to self- 
consistent values of the gap parameters at 8 K of Ai = 
1.25 meV and A3 = -0.53 meV from Eq.ElandEl These 
parameter values yield Tc= 15 K, which is in reasonable 
accord with the experimental value of 26 K. Note that 



a leading edge gap is clearly seen near the hotspot in 
(b). This gap is a direct consequence of the broadening 
of the spectrum due to small angle scattering and the as- 
sociated residual spectral weight at the Ep, even though 
this spectral weight is quite small due to the presence of 
a sizable pseudogap in the spectrum. 

Fig. 4(a) shows variations in the leading edge shift 
Ashift and its non-monotonic and anisotropic nature more 
clearly. Here the minimum computed value of Ashift (blue 
line) is plotted along various directions given by the an- 
gle (/), where (j) — Q refers to the antinodal direction and 
(j) — 45° to the nodal direction. The actual momentum 
points at which the computed values of Aghift are plot- 
ted, shown in the inset, lie close to the non-interacting 
FS. The corresponding experimental values (red dots) are 
in good accord with the theoretical results. In particu- 
lar, the gap reaches a maximum value of about 2.5 meV 
around (p = 18° near the hot-spot in both theory and 
experiment. The monotonic d-wave gap obtained by set- 
ting A3 = in the calculations is also shown for reference 
(dashed line) to highlight the non-monotonic gap varia- 
tions in the present case. 

Our explanation of nonmonotonic gap variations 
should be distinguished sharply from the notion that one 
simply observes two different gaps in the experiments, 
i.e. an SC gap along the pockets and an AFM gap near 
the hotspots. This alternative picture suffers from the 
problem that the AFM gap in the normal state of the 
cuprates is far too large, and moreover, it is not clear 
how the leading edge gap comes about since there are no 
states at the Ep away from the region of the two small 
FS pockets. Interestingly, the two band model study of 
Ref. llSl would require an anomalously small value of t to 
fit the experimental gap values in PLCCO. 

Third harmonic gaps have been found not only in 
electron-doped PLCCO0 and NCCOQ, but also in 
hole-doped (underdoped) Bi2Sr2Ca„_iCu„02n+4 {n = 
1-3) (BISCO)!!^,!!^. Interestingly, the third harmonic 
contribution in electron-doped cuprates possesses an op- 
posite sign and its size is substantially larger than in 
the hole-doped case, leading to the nonmonotonic gap 
variations discussed above. Our study provides insight 
into this behavior. Separate electron and hole pockets in 
PLCCO and NCCO are suggestive of two-band|22J su- 
perconductivity as in the case of MgB2 where interband 
pairing greatly enhances the Tc|23]. However, the form 
of our Hamiltonian in Eq. ^ does not involve a sepa- 
rate interband pairing gap. Fig. Efb) shows that A^j: 
possesses the correct symmetry and that its size is maxi- 
mal in the interval between the two pockets just like the 
leading edge gap Aghift- With this in mind, we propose 
that the third-harmonic gap A,^ can be looked upon as 
playing the role of interband pairing gap in PLCCO and 
NCCO. In fact, the 'hot spot' between the two pockets is 
associated with strong AFM fluctuations in a quasi-two- 
dimensional svstem|4l lll|. which in turn are responsible 
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FIG. 4: (color online) (a) Computed leading-edge gap (blue 
line) Ashift at 8 K as a function of the FS angle <p is com- 
pared with the corresponding experimental results (red dots) 
on PLCCO 3l- A-D are momentum points discussed previ- 
ously in connection with Fig. 3. Green line gives the mono- 
tonic gap variation at 8 K in the absence of the third harmonic 
term (i.e. A3 = 0). The inset shows (open circles) the actual 
momentum points for various directions <p where the com- 
puted gap values are plotted, (b) Form of the third harmonic 
contribution A^g (magneta curve) to the gap (see Eq. 3) is 
shown schematically on top of the FS map of Fig. 1(b). 



for coupling the UHB and the LHB, and thus coupling 
the phase of the order parameter on the two FS segments. 

The fact that the nonmonotonic gap in our model 
is intimately associated with the AFM pseudogap in- 
dicates that superconductivity arises in an AFM back- 
ground. These results support the picture of electron 
doped cuprates being uniformly doped AFM metals, even 
in the superconducting state, rather than being phase 
separated AFM and SC domains. Note however that 
when fluctuations are added, the present AFM gap be- 
comes a short-range order pseudogap, with a possible 
Neel ordering transition at lower temperatures due to in- 
terlayer coupling. Open questions remain as to whether 
superconductivity coexists with this residual 3D Neel or- 
der, or whether the onset of superconductivity destroys 
either the 3D order, or even the 2D (T=OK) long-range 
AFM order. 
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